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SUMMARY 


The theory of small disturbances, -i s applied to the calculation 
of the pressure distribution and drag of a closed body of revolution 
traveling at supersonic speeds. It is shown that toward the rear 
of the body the shape of the pressure distribution is similar to 
that for subsonic flow. For fineness ratios between 10 ahd lf> the 
theoretical wave drag is of the same order as probable values of 
the frictional drag . ... 


INTRODUCTION 


Methods for calculation of the flow over a body of revolution 
traveling at supersonic velocities have been known for some time. 
(See references 1 and 2.) Investigations along these lines have, 
however, been confined chiefly to bodies having the form of 
artillery projectiles. Such bodies, because of their blunt forms, 
show relatively high drags and are thus not suited for use on 
high-speed aircraft. The drag of slender bodies and the effects 
of fairing the rear of these bodies are therefore of considerable 
interest in connection with the problem of flight at speeds above 
the speed of sound. 

In view of the interest in possible aeronautical applications 

it was thought worth while to apply the known methods to a 

particular case of a closed body having both a tapered nose and a 

tapered tail. Slender shapes described by .the potation of 

parabolic arcs were chosen and the resulting pressure distributions 

over the surfaces as well as along the axes behind the bodies were 

.calculated. The results are compare^ with those obtained for 

similar shapes in an incompressible fluid -and el so in one case 

with a two -dimensional body having a similar cross section. 

* 
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SYMBOLS 
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V 

a 

M 

B’ » 
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p , r 

q . 

Ap 

0 

l 

R 

,d 

L 

s max 

D 

C D* 


z Cartesian coordinates 

undisturbed fluid velocity- 
speed of sound in fluid 
Mach number (V/a) 

\( M 2 - 1 


T 


Cp ' 

"■tn 


VOl 


density of fluid 
dynamic pressure ^ 
pressure increment 

* * 

velocity potential of single source 


1 


velocity potential of continuous distribution 
of sources along x-axis 

•abscissa of individual sources’. 

radius of body 

maximum diameter of body 

length of body 

maximum cross-sectional area 
drag 

drag coefficient based on maximum cross-sectional 
area (D/qSj^) 

drag coefficient based on: 2/3 power of 

volume f- — 2 J 2 

t q( Volume) ' 


c, C 


constants 


V 
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METHODS OF CALCULATION 

I 


The method used herein follows closely that of reference 1. 
Figure 1 shows the shape of the body and the orientation of the 
axes. The disturbance produced by the body is assumed to be small 
and the flow isentropic so that the linearized equation for the 
potential of the disturbance velocities • 0 will apply. This 
equation is (see reference 2) = 

C 1 - « e ) '-.fa '+ ?yy + - 0 (1) 

As in the case of an incompressible fluid the flow over 
the body can be obtained by the addition of flows due to an 
infinite number of sources distributed along the axis. The 
potential of a single source in a supersonic stream is 


where 


and 


Figure 2 shows 

source compared with those for a source in an incompressible 
flow. In the case of a source in an incompressible flow the 
equipotential surfaces are spheres, given by the expression 


CV 


+ z 2 ") 

-C.V- ■ -- 

\/* 2 - B 2 r 2 . 


( 2 ) 


B - \/ M 2 - 1 


r ■ \/ y 2 +' z 2 

the equipotential lines for the supersonic 


i*o- 


SlL. 




T~ 2~ 2 

+ y + z 


( 3 ) 


In the supersonic case the equipotential surfaces are hyperboloids 
of two sheets contained v/ithin the Mach cones. Although the 
mathematical expression has values in two cones, one ahead of and 
one behind the source, only the values beliind have physical 
significance . 
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It mil be noted that the distribution of velocities along 
the x-axis is the same for the supersonic source (equation (2)) 
as for the subsonic or incompressible source (equation (3)). 

Since -the forward cone is to be disregarded in the supersonic case, 
however, it .is found that the coefficient G in equation (2) 
must be doubled in order to produce the same flux, or intensity, 
as equation (3) The result is that the velocities along the 
axis behind, a; supersonic, source are exactly twice those of a 
subsonic source having tho same intensity.- • 


The sources and sinks are assumed to be continuously 
distributed with intensity ' 27 tVf(£) • per unit length along the 
x-axis from -1 to 1. The abscissas of the individual sources are 
denoted by 5,. Positive values of f(£) denote sources and 
negative values denote sinks.- By adding the potentials due to 
tho single elementary sources f(2,)<S,. the resulant flow 


0 




JT(.£).4S_ 

£) 2 - B 2 r 2 

... 


(M 


is obtained. 


The problem is to determine a source distribution in such a 
way that — 



( 5 ) 


where dR/dx gives the shape of the meridian curve of the body 
of revolution. It is shown in reference 1 that to a first 
approximation for a slender body the source strength' is pro- 
portional to the rate of change of . the cross section of the body, 
that is 


fU)- K g .... t6) 

A similar approximation can be app r licd to' obtain the source 
distribution for a body in subsonic flow. The distributions are, 
in fact, 1 , the same in the two cases' with the ; exception that in the 
supersonic flow the value of f(x) must be doubled to account 
for the elimination, of the flux through the forward cones. 
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By choosing f(x) = 
following expression for 


E 


c(x^ - x) and solving equation (6), the 
R was obtained: .,G^' 




' ( T ) 




This expression may “be recognized as the equation of a surface 
obtained by revolving a parabolic arc about its chord. The 
fineness ratio of the body is determined by the value assigned 
the factor c. 

On substituting c(<S^ - &) for f(£) in equation (*!•)* the 
velocity increment 30/ox at point (x>r ) is found to bo (see 
equation (9*5), P*39> reference l) 


7 x-Br 

- v/ Zl&l 


-1 \/(x - i f - B 2 r 2 


-Br 


■ cV 
V 


: 13S 1-hJ.l 

^ - 5) 2 - 




= cV 


; (3 r 9*) . \/{ x + l) 2 - B 2 r 2 ' 


x+l> 


+ ( 3x 2 - 1 + | B 2 r 2 ) cosh ( Tr 


over the body. 

Along the axis behind the body the integration gives 


cY 


(3x 2 - 1) los'arHr " 6x 


(S) 


(9) 


The prossur*. coefficients wore calculated by tho .-formula 


^=4- g 

q V u 


( 10 ) 
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RESULTS 


Qalculations have been made for a Mach number of I.I 4 and for 
three thickness ratios d/L of 0,066?,, 0,10* and corresponding 

to fineness ratios .I/d of l£, 10, and 6, 67 , respectively. The 
results' for the thre.4 bodies are shown in figures 3, h , and ir- 
respectively, and are compared with the theoretical pressure-— ■ 
distributions over these bodies in an incompressible, fluid. A 
discussion of the errors involved in the linear theory and the 
variation of the pressures with Mach number will- be found in 
ref er4rice ' 3 • , . .... 

Comparison of the distributions in a compressible 'fluid and 
in an incompressible fluid shows a certain similarity, especially 
toward the rear ef the body. The effect of supersonic speed 
appears to be similar to the effect of .a lag inasmuch as the _ 
negative pressure peak and the region of pressure recovery are 
displaced rearward. The pressures alone the axis behind the bod y 
a re .just twice those produced bv an incompressible fluid, as may - 
be seen bv referring to the velocity field of a~ single source ^ 

The results obtained herein for the three-dimensional body 
are in marked contrast to the results obtained for two-dimens .1 onal 
bodies, or wing sections, having similar cross sections 1 .- As is 
we ll known, in the two-dimensional case no pressure recovery 
t akesolace at supersonic speeds, th e pressure at"~a poihlT'belng 
det ermined so leivlhv -the .innl i nation of the surface -at - that p oint 
so that positive pressures occur wherever the cross section is 
expanding and negative pressures occur wherever the cross 
section is diminishing . Figure ,6 shows tho comparison of the 
te/o-dlmens ibnal and three-dimensional bodies for the , 0*10 
thickness ratio,' 

The essential difference between the two- and throe- , ... 
dimensional flows corresponds to the difference noted by Lamb 
(reference ii) between the characteristics of a plane sound wave 
and an axially symmetrical wave diverging from a center. As 
noted by Lamb, the piano wave, which corresponds in the present 
case to the flow produced by the two-dimensional wing section, 
is propagated indefinitely without change of formj whereas tho 
axially symmetrical wave, which approximates that produced by an 
element of the slondor body of revolution, does not follow the 
form of- the disturbing motion but loaves a "tail" of diminishing 
intensity and indefinite extent. Thus tho wing soot.i on leave s 
no pressure disturbance in its wake, whereas tho axially 
symme trical body is followed by an indefinite region of positiv e 
pressure. Integration or the axial components of the prcs'sui*6s 
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acting on the tody, however, shows that the positive pressure at 
the rear of the three-dimensional body is sufficient to cancel, 
only a small fraction of the total pressure or wave drag. 

The wave-drag coefficients based on the maximum frontal area 
were found to be O.OI 4 . 9 , 0.11, and 0.2U for. the. bodies with thickness 
ratios of 0. 0667, _ 0,10, and ©.Irrespectively. - Figure 7 shows 
a comparison of these values with the wave drags of corresponding 
two-dimensional wing sections. It will be noted that the wave 
drag of the fuselage form is approximately proportional to the 
square of the thickness ratio. - . .. 

An approximate estimate of the total drag of a body may be 
obtained by adding values of the frictional drag to the wave drag. 

A rather complete treatment of the frictional drag of bodies of 
revolution at subsonic speeds is available from reference 5. By 
use of values from reference 5 corresponding to a fully turbulent 
boundary layer and a Reynolds number of 10®, the following estimates 
of the total drags of the bodies were obtained: 


. * ' . • v 

Thickness ratio 


Cd vo1 

0.0667 

0 . 1 U 

0.032 

.10 

. 17 .... 

■ .051 • 

.15. 

.29 

.11 


where Cj-j is the drag, coefficient based on the frontal area and 

On is the drag coefficient based on the 'volume of the body to 
vol 

the 2/3 power. , ' • ‘ 

The drag of a given volume is -an .important criterion in the 
case of an airplane fuselage and it will be of interest to compare 
these values with a typical value attainable at subsonic speeds. 
For a Reynolds number of 10° and a turbulent boundary layer, the 
best value given by Young (reference 5) corresponds to a thickness 
ratio of 0.2 and is approximately 

CL = 0.016 
D vol 

CONCLUDING REMARKS 

The theoretical pressure distribution over a closed body of 
revolution traveling at supersonic velocities shows a pressure 
recovery at the rear of the body similar to that occurring at 

7 
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subsonic speeds . The extent of the region of 'positive pressure 
is, however, not sufficient to have a pronounced effect on the 
wave drag. It appears to be necessary to use extremely slender 
shapes to obtain total drag values comparable to those of a 
conventional airplane fuselage at subsonic speeds. For fineness 
ratios between 10 and If? the theoretical wave' drag is' of the seine 
order as probable values of the frictional drag. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., July 8 , 19k& 
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Figure 1.- Orientation of axes and shape of body. 
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Figure 2.- Equipotential lines for a supersonic source compared with those 

for a source in an incompressible flow. 
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Figure 3.- Calculated pressure distributions for a thickness ratio 

of 0.0667. 
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Figure 4 .- Calculated pressure distributions for a thickness ratio 

of 0.10. 
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Figure 5.- Calculated pressure distributions for a thickness ratio 

of 0.15. 
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